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1 .  Report  Susrarv 


Deep  atomic  hole  states  constitute  a  regime  that  has  only 
recently  Decome  amenable  to  detailed  investigation,  thanks  to 
large-scale  computers  and  such  energetic  probes  as  svncnrotron 
radiation  and  accelerated  ion  beams.  The  physics  of  atomic  inner 
shells  is  peculiar  in  several  respects:  ( i )  The  energies  are 
high;  e.g.,  the  Is  binding  energy  reacnes  100  keV  at  Z=87.  Conse¬ 
quently,  strong  relativistic  ana  quantum  electrodynamic  effects 
arise,  (ii)  Transitions  are  mostly  radiationless,  proceed  through 
many  channels,  and  are  often  very  strong;  hole-state  lifetimes 
thus  become  as  short  as  attoseconds,  with  widths  measured  in  tens 
of  electron  volts.  The  dynamics  of  inner-shell  transitions  conse¬ 
quently  exhibit  unusual  features. 

The  central  purpose  of  our  research  program  is  to  explore 
tnese  peculiar  properties  of  energetic  inner-shell  processes.  We 
have,  in  the  past,  paid  particular  attention  to  relativistic  and 
quantum-electrodynamic  effects  and  have  performed  large-scale 
state-of-the-art  computions  of  hole-state  energies  and  transition 
rates.  In  the  course  of  this  work,  we  have  developed  the  first 
general  computer  code  for  the  relativistic  calculation  of  Auger 
rates . * 

At  the  present  time,  our  primary  aim  is  to  gain  a  better 
understanding  of  the  dynamics  of  inner-shell  processes,  particu- 
*  *ariv  under  threshold  excitation.  It  has  been  known  for  some 

years  that  here  the  traditional  two-step  model  breaks  down,  i.e., 

I  one  cannot  treat  the  decay  of  a  deep  inner-shell  hole  state  as 

separate  from  its  production.  The  concept  breaks  down  that  there 
is  a  relaxation  pnase  during  which  the  electrons  in  the  excited 
atom  adjust  themselves  to  the  new  ->otential  and  "forget"  how  the 
state  was  produced.  Instead,  threshold  excitation  and  decay  take 
place  in  a  single  second-order  transition,  which  in  essence  can 
'  be  thought  of  as  the  inelpstic  analog  of  resonance  fluorescence. 

^  A  few  years  ago,  we  were  aoie  to  demonstrate  the  single-step 

nature  of  inner-shell  excitation  and  radiationless  deexcitation, 
dubbed  the  "resonant  Auger  Raman  effect."2  This  effect  results  in 
-J  resonant  emission  of  sub-natural  linewidth  radiation. 


r 

c 

r 


| 

/ 

/ 

V 

■ 


With  the  single-step  nature  of  threshold  excitation-deexci¬ 
tation  established,  and  the  two-step  model  long  proved  valid  in 
the  asymptotic  limit,  the  pressing  fundamental  question  is:  how 
does  nature  link  these  seemingly  contradictory  regimes?  We  have 
devoted  several  years  to  an  investigation  of  this  question,  and 
to  a  search  for  a  unified  description  that  would  encompass  the 
entire  range  of  excitation  energies. 

Related  questions  with  which  we  have  been  concerned  are  (i! 
the  extent  and  nature  of  many-body  processes  associated  witn 
atomic  inner-shell  excitation,  (ii)  a  correct  relativistic  des¬ 
cription  of  inner-shell  ionization  by  slow  charged  particles,  and 
of  subsequent  x-ray  emission,  and  (iii)  a  search  for  an  expiana- 


tion  why  theory  almost  Invariably  predicts  much  too  short  a 
lifetime  for  certain  atomic  hole  states,  such  as  [2s]  . 

j 

We  have  pursued  an  integrated  theoretical  and  experimental 
approach.  The  theoretical  aspect  has  a  strong  numerical  compo¬ 
nent.  Major  computations  are  being  performed  at  the  NASA  Ames  i 

Research  Center,  whicn  has  a  CRAY  1  computer.  A  professional 
programmer  works  half-time  for  us  at  NASA  ARC.  The  experimental 
work  is  primarily  conducted  in  the  Stanford  Synchrotron  Radiation 
Laboratory  (SSRL),  where  we  operate  an  on-line  gas-pnase  electron 
spectrometry  system.  Data  are  transferred  by  magnetic  tape  to 
Oregon  for  analysis  on  campus.  j 

The  major  result  of  the  year  is  the  successful  outcome  of 
our  investigation  of  the  dynamics  of  inner-shell  photoionization 
and  radiationless  deexcitation  (see  Sec.  2.1).  Detailed  measure¬ 
ments  of  the  Auger  resonant  Raman  effect  in  the  L3  subshell  of 
xenon  and  of  post-collision  interaction  (PCI)  shifts  of  the  Lg- 
m4m5  '  G4)  Auger  line  were  completed  and  successfulv  explained  in 
the  framework  of  resonant  scattering  theory.  For  this  purpose, 
the  first  fui^y  quantum-mechanical  theory  of  post-collision  in¬ 
teraction  was  constructed.  The  result  is  a  consistent,  unified 
description  of  excitation/deexcitation  of  deep  hole  states,  that 
incorporates  the  resonant  Raman  effect  at  threshold  and  asympto¬ 
tically  turns  into  tne  two-step  model,  with  PCI  providing* the  ! 

link  between  the  two  extremes.  This  result  culminates  five  years 
of  intensive  work.  Crucial  contributors  were  Teijo  Aberg  of  the 
Helsinki  University  of  Technology,  who  spent  the  fall  quarter  in 
our  group  as  a  visiting  professor  and  developed  the  resonant 
scattering  formulation,  and  Mau  Hsiung  Chen  of  the  Lawrence 
Livermore  Laboratory,  who  solved  the  extremely  difficult  two-  l 

electron  continuum  wave-function  problem. 

A  related  important  result  of  this  year's  work  concerns  a 
many-body  effect  in  inner-shell  ionization.  A  synchrotron-radia¬ 
tion  study  of  double  photoionization,  through  measurements  of 
Auger  satellites,  permitted  us  to  demonstrate  for  the  first  time  I 

the  theoretically  predicted  difference  between  the  dependence  of 
shakeup  and  shakeoff  probabilities  on  the  ohoton  energy  near 
threshold  (cf.  Sec.  2. 2). 3 

We  have  utilized  the  large  computing  power  of  the  NASA  ARC 
facility  in  order  to  construct  tables  of  K-  and  L-shell  . 

crossections  for  ionization  under  impact  by  slow  protons,  on  the 
basis  of  a  fully  relativistic  theory  that’we  develODed 
previously4  (Sec.  2.3).  *  } 


Progress  was  achieved  in  a  protracted  theoretical  study 
designed  to  understand  why  calculated  lifetimes  of  certain  inner- 
shell  hole  states  (e.g.,  in  the  2s  subshell)  tend  to  be  very  much 
shorter  (by  a  factor  of  2  or  3)  than  measured.  This  problem  is 
perplexing  because  calculations  of  other  lifetimes 
tend  to  come  out  well;  it  is  important  because  the 
such  calculations  enter  into  the  design  of  schemes 
wavelength  lasers.  Added  to  previous  studies  of  the 


(or  widths) 
results  of 
for  short- 
effects  of 
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relaxation  and  intermediate  coupling,5,6  we  have  now  calculated 
the  effects  of  final-state  channel  mixing.  The  interchannel  in¬ 
teraction  is  found  to  have  a  major  influence  on  transition  rates 
(Sec.  2.4). 

A  calculation  of  x-ray  emission  rates  for  transitions  to 
atomic  M-shell  vacancies  has  been  completed,  with  particular 
emphasis  on  the  results  of  the  choice  of  gauge.'  This  is  the 
first  sytematic  relativistic  calculation  of  such  rates  employing 
a  nonlocal  potential  (Sec.  2.5). 

Our  results  on  dynamical  and  many-bodv  aspects  of  inner- 
shell  transitions  urgently  call  for  extension  of  these  studies. 
We  plan  to  carry  out  further  calculations  of  post-collision 
interaction  to  encompass  other  experimental  data  obtained  re¬ 
cently,  and  to  perform  new  high-resolution  measurements  in  the 
threshold  region  where  additional  data  can  provide  critical  tests 
of  the  resonant  scattering  theory.  For  this  purpose,  we  hope  to 
gain  access  to  the  new  54-pole  wiggler  line  at  SSRL,  which  pro¬ 
vides  substantially  enhanced  x-ray  flux.  Theoretical  results  on 
inner-shell  transition  rates  and  ionization  cross  sections  also 
invite  significant  extensions,  building  on  the  recently  attained 
insights. 
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2 .  Research  Results 


2 . 1  Dynamics  of  atomic  inner-shell  processes 

In  threshold  excitation  of  short-lived,  deep  atomic  inner- 
shell  hole  states,  ionization  and  decay  cannot  be  treated  as 
distinct  processes.  Directly  at  threshold,  photoionization  and 
radiationless  deexcitation  occur  in  a  single  process,  the  reso¬ 
nant  Raman  effect.1  Above  threshold,  excitation  and  deexcitation 
are  still  linked  by  post-collision  interaction  (PCI),  in  which 
the  Auger  decay  takes  place  under  the  influence  of  the  Coulomb 
field  of  the  receding  photoelectron  and  some  of  the  photoeiec- 
i  tron's  energy  is  transferred  to  the  Auger  electron.  Measurements 

|  of  the  Auger  spectrum  excited  near  threshold  can  therefore  probe 


' 

i 


the  dynamics  of  inner-shell  transitions  and  provide  a  stringent 
test  of  the  semic lassical  PCI  theory2  whicn  has  been  shown  to  be 
valid  in  the  limit  of  long  hole-state  lifetimes.3  The  first  fully 
quantum-mechanical  calculation  of  PCI  becomes  possible,  based  on 
resonant  scattering  theory,*'5  since  the  intermediate-state  sum¬ 
mations  can  be  restricted  to  states  associated  with  the  initial 
inner-shell  hole. 

We  have  performed  an  investigation  in  which  highly  monochro- 
matized,  hard  synchrotron  radiation  was  tuned  through  the  L3 
threshold  of  Xe  and  the  L3-M4M5  (1G4)  Auger  decay  including  the 
nd  (n>5)  spectator-electron  satellite  distribution  was  observed. 
The  peak  positions  and  widths  of  the  n=5  and  n*6  spectator  lines 
were  measured  as  functions  of  photon  energy.  The  evolution  of  the 
M4M3  double  photoionization  peak  into  the  Auger  diagram  line  was 
found  to  be  strongly  blended  with  the  n>6  spectator  satellites  in 
the  excitation-energy  range  from  Eexc=-0.5  eV  to  +5  eV  with 
reference  to  the  ionization  threshold,  producing  a  large  apparent 
PCI  shift  (-3  eV).  We  compare  this  result  with  predictions  from 
the  resonant  scattering  theory  at  Eexc*+3  eV.  Above  Eexc*+5  eV, 
the  PCI  shift  is  solely  due  to  distortion  of  the  diacr'am- 1  ine 
snape,  so  that  comparison  with  both  the  resonant  scattering  and 
semiciassical  theories  becomes  possible. 

In  the  experiment,  performed  in  the  Stanford  Synchrotron 
Radiation  Laboratory,  x  rays  from  an  8-pole  wiggler  operating  at 
14  kG  were  focussed  onto  a  Xe  jet  by  a  Pt-coated  doubly-curved 
toroidal  mirror.  A  narrow  energy  band  (-0.6  eV)  was  selected  with 
a  tunable  Si  (111)  (1,-1)  double-crystal  monochromator.  The  ener¬ 
gy  spectrum  of  electrons  emitted  from  the  Xe  target  was  measured 
with  a  double-pass  cyl indr ical -mirror  analyzer.1  We  take  the 
asymptotic  Auger  diagram-line  energy  6A°=3 , 365 .9+0.2  eV  as  refe¬ 
rence  for  the  PCI  shift.  Typical  Auger  spectra  excited  at  Eexc=-2 
eV  and  +1.S  eV  are  shown  in  Fig.  1. 

In  the  resonant  scattering  theory*'5  of  the  PCI  effect,  a 
consistent  treatment  of  threshold  phenomena  is  achieved  by  inter¬ 
preting  the  Auger  transition  [ni 1 ^ ] * [nf 1 f , n<= .  1  f ,  ]  as  a  resonance 
in  the  double  electron  pho to ionizat ion  of’ the  nf^f  and  nf,lf, 
subshells.6  As  the  incident-photon  energy  approaches  the  ioniza¬ 
tion  energy  IA  of  the  i^l^  threshold,  spectator  transitions  of 
the  type  [n^ljjnl*  [nf  1^  ,nf  .1  *,]nl  start  to  appear  owing  to  n.Ii? 
nl  excitations.  At  the  threshold,  the  spectator  lines  blend  into 
the  rapidly  increasing  double-ionization  peak,  for  which  the 
cross  section  at  large  excitation  energies  Eexc=<0-Ij.  becomes 
the  product  of  the  n^l^>£l  single-electron  photc ionization  cross 
section  and  the  Auger*  decay  probability  PC**®)  for  the  (n^jj 
> [n^l f , n^r  1  f.  ] 2S+1Lj  transitions  with  energy  £A  .  Close  to  thres¬ 
hold,  the  Auger  decay  cannot  be  treated  as  an  independent  process 
because  the  emission  of  the  Auger  electron  is  affected  by  transi- 
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FIG.  1.  Xenon  L3-M4M5  (1G4) 
Auger-electron  spectrum,  excited  with 
synchrotron-radiation  photons  of 
energy  2.0  eV  below  the  Xe  L3  binding 
energy  (top)  and  1.5  eV  above  the  Xe 
L3  binding  energy  (bottom).  Feature  A 
consists  primarily  of  nd  (n>6)  satel¬ 
lites,  with  a  slight  admixture  of 
the  PCI-shifted  diagram  line;  feature 
B  is  the  5d  spectator  satellite,  and 
feature  C  consists  of  the  PCI-shifted 
diagram  line,  with  an  admixture  of  nd 
(n>6)  satellites. 


tions  from  the  single-  to  double-hole  n(£)l  states.  The  double 
photoionization  cross  section  is 


‘I  tracts. 


( i ) 


where  N(u-u>»)  is  the  distribution  function  of  the  incident  pho¬ 
tons  (normalized  to  unity),  dA  is  the  Auger-electron  energy,  and 
I**i  is  the 


‘ft' 

have 


ionization  energy "of  the  final  double-hole  state.  We 
placed  the  slowly  varying  factor  u>oOrgA)  outside  the  integ¬ 
ral.  If  the  continuum  wave  function  (t>  is  replaced  by  the  dis¬ 
crete-state  wave  function  ln>  in  the  overlap  matrix  element 
<<  !t'>,  the  double-ionization  cross  section  reduces  to  the  cross 
section  d(TL/d&A  for  emitting  an  Auger  electron  and  placing  a 
spectator  electron  in  the  excited  state  nl.  In  Eq.  (1),  we  have 

I  ^><tl|r||0>ir 
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Eln><n  Hr  (|  0>  C 
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(2) 


where  each  |t(n')>  must  be  evaluated  in  the  field  of  the  singly 
ionized  atom  with  an  n^l^  hole,  in  contrast  to  the  final  |£(n)> 
double-hole  states.  In  Eq.  (2),  T.  is  the  total  decay  width  of 
the  [n^lj]  hole  state,  and  <  Z  (n1)  ||  r||  0>  is  the  reduced  radial 
dipole  matrix  element. 


as  a 


Close 

function 


to  threshold,  the  Auger-electron  emission 
of*A  is  a  measure  of  the  cross  section 


intensity 


(3) 


which  must  be  convoluted  with  the  spectrometer  window  function 
for  comparison  with  the  experimental  data.  If  the  final  double- 
hole-state  lifetime  Pf-1  is  taken  into  account,  the  Dirac  delta 
function  in  Eq.  (1)  is  replaced  by  the  normalized  density  fun¬ 
ction  l  /1JT>C(  <«*- 6  -  **  -  Iff'  )v  ♦  It  follows  from  Eqs. 

(1)  and  (2)  that  for  £^,.>0  the  discrete  part  of  the  cross  sec¬ 
tion  (3)  rapidly  vanishes  compared  with  the  continuum  part  dr^./d*.. 

In  the  n=>0  limit,  the  corresponding  overlap  integral 
<£\V>  can  be  approximated  by  taking  the  overlap  between  two  WKB 
wave  functions,  corresponding  to  ■£  k  ^  -Eexc  + i  (<”\  /  2  )  ♦  V  .  ( R )  and 
3kf 2=Eexc-<A+ £A°  +  Vf(R)  ,  respectively.  we  have  Vf  =  2Vi=R“*,  where 
the  angular-momentum  exchange  between  the  photoelectron  and  the 
Auger  electron  is  neglected.  All  semic laosical  models  of  PCI  seem 
to  be  variants  of  the  WKB  approximation,  limited  to  the  verv- 
high-n  and  continuous-^  region.  Qualitatively,  both  Eqs.  ( 1 )  —  ( 3 ) 
and  the  semi c 1 assi ca 1  approach  lead  to  a  transfer  of  intensity 
from  the  low-energy  to  tne  high-energy  flank  of  the  Auger-elec¬ 
tron  peak,  compared  with  a  Lorentzian  shape,  resulting  in  a 
positive  PCI  shift  .  However,  whereas  in  the  semiclassical 
approach  the  £  »  -tA°5Eexc  regime  is  approached  by  analytic 

continuation,  tne  quantum-mechanical  theory  leads  to  a  spectator¬ 
line  structure,  important  in  practice  for  Eexc?2^i- 

The  derivation  of  Eq.  (1)  from  the  general  transition 
matrix  formula4  is  based  on  a  number  of  approximations.  The  most 
severe  of  these  is  probably  the  factorization  of  the  many-elec- 
tron  Hamiltonian  matrix  element  that  involves  the  final  two- 
electron  scattering  wave  function  into  the  one-electron  overlap 
element  <1  If >  and  the  Auger-electron  probability  amplitude. 

^  The  cross  section  (3)  was  calculated  for  the  L3-M4M5 

(1G4)  transition  with  the  use  of  Hartree-Fock  (HF)  wave  fun¬ 
ctions.  Calculation  of  the  continuous  part  of  the  function  (2) 
required,  for  each  Eexc  (3,  10,  and  20  eV)  ,  2,000  continuum  wave 
functions  l£d>  in  steps  of  0.1  eV.  These  wave  functions,  as  well 
as  the  final -state  continuum  wave  functions  IX.d>,  were  generated 
in  the  configuration-average  frozen  HF  core,  with  the  exchange 
interaction  taken  into  account.  The  summation  over  the  discrete 
intermediate  states  was  found  to  be  negligible  in  all  three 
cases . 

The  -2.3-eV  width  of  the  spectrometer  window  allows  us  to 
resolve  only  the  5d  spectator  line  in  the  range  -8e V<Eexc<+5e V , 
at  5d) =Eexc+ 10 . 5  eV,  from  the  rest  of  the  structure  which 
appears  as  a  single  peak  at  £<A(5d).  Up  to  Eexc=0  eV,  A  ap¬ 
pears  to  follow  a  linear  dispersion  law,  A=Eexc+4.0  eV,  which 
indicates  that  it  is  mostly  6d.  The  linear  dependence  of  A(5d) 
and  £(6d)  on  Eexc  is  plotted  in  Fig.  2;  the  resonance  energies 
are  indicated  at  which  the  spectator  lines  gain  maximum  intensi¬ 
ty.  According  to  Eqs.  ( 1 )  —  ( 3 )  ,  these  resonances  occur  at  = 

C  Aw 
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or  £(nd)*£n-tn.  The  measured  energies  Cn  agree  well  with  re¬ 
laxed  HF  calculations  for  n=5  and  6. 
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FIG.  2.  Xenon  L3-M4M5  (*G4)  Auger-electron-peak  energies, 
as  functions  of  excitation  energy.  Experimental  data  are  indi¬ 
cated  by  solid  circles.  Diamonds  represent  the  quantum-mechanical 
shift  &  of  the  diagram  line,  while  the  solid  curve  indicates  the 
semiclassical  shift. 


The  measured  width  of  the  5d  and  6d  spectator  lines  is 
4. 0+0. 5  eV,  independently  of  Eexc.  This  is  clearly  narrower  than 
the  natural  width  of  the  L3-M4M5  (1G4)  line  ( f^+Hj  *4 . 2  eV),7 
convoluted  with  the  2.3-eV  wide  spectrometer  window,  i.e.,  5.3 
eV.  This  resonance  narrowing  is  in  accordance  with  Eqs.  <  1 )  —  ( 3  ) 
which,  in  analogy  with  the  x-ray  resonant  Raman  case,5  predict 
<Afl“n/d£A  oC  where  N(n-«0):s  0.6  eV  wide  in  our  experiment. 

With  the  sum  rule  J | < £ 1 1 ' >  |  2d£ = < c'  I  C  >  ,  where  summation 
over  discrete  states  is  included,  it  is  possible  to  estimate  how 
much  of  the  low-energy  structure  represented  by  peak  C  in  Fig.  1 
is  due  to  spectator  states.  The  result,  in  agreement  with  the 
gradual,  slow  disappearance  of  the  5d  spectator  line,  is  that  for 
2exc?2*"i'  the  contribution  from  the  discrete  final  states  cannot 
be*  neglected.  Hence  Eexc=+3  eV  represents  an  intermediate  case, 
in  which  there  is  an  apparent  shift  due  to  the  spectators.  Expe¬ 
rimentally  we  find  A=3.0+0.5  eV,  whereas  theory  predicts  &=2.8  eV 
after  convoluting  with  the  final -state  density  and  the  spectrome¬ 
ter  window  function. 
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In  the  calculation  of  d0-n/dd.A,  the  discrete  portion  of 
complete-set  HF  computations  was  approximated  by  quantum  defect 
theory  for  both  n'  (intermediate)  and  n  (final)  states  >8.  The 
summation  over  intermediate  continuum  states  was  found  to  be 
important,  indicating  that  there  is  appreciable  recapture  of  the 
photoeiectron  into  Rydberg  states. 


For 


Eexc55 
diagram  line 


eV,  the  quantum-mechanical  shape  of  the  Auger 
agrees  closely  with  that  obtained  from  the  semiclas- 
sicai  Niehaus  theory2  through  numerical  integration.  After  convo¬ 
lving  the  theoretical  shapes,  we  find  the  shifts  ^theorv*2,0 

^expt  =  2,5-0,5  eV  for  Eexc* 10  eV-  and  eV 

1 . 9  +  0 . 7  eV  for  E  -20  eV.  Si 

the'”'”  * 


eV 

a. 


caff 


A*i 

y  why  the 


1  eV'  and  Atheorv‘ 
t  remains  to  be  shown"  analyti- 

Niehaus~£fieory, *  if  the  stationary-phase  approxima¬ 
tion  is  not  invoked,  leads  to  almost  Identically  the  same  re¬ 
sults  as  the  quantum-mechanical  approach.  In  Fig.  2,  the  solid 
curve  for  the  diagram-line  energy  represents  the  function 
£=&<Eexc),  obtained  from  the  Niehaus  theory  and  convoluted. 


Similar  measurements  as  here  described  have  been  perfor¬ 
med  on  other  wide  inner-shell  hole  states.^  Computations  are  in 
progress  to  apply  the  resonant  scattering  theory4'5  to  these 
cases  as  well.  A  preliminary  report  on  this  work  has  been 
submitted  for  publication  (cf.  Sec.  3,  Item  12  of  this  Report). 
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2.2  Threshold  double  oho to ionization  of  argon  with  synchrotron 
radiation. 

In  atomic  inner-shell  photoionization,  multiple  excitation 
processes  occur  with  significant  probability.  The  resulting  final 
states  are  approximately  described  by  configurations  formed  by 
removal  of  a  core  electron  and  excitation  of  additional  electrons 
to  higher  bound  states  (shakeup)  or  to  the  continuum 
(shakeof f ) . 1-3  Such  multiple  excitation  Drocesses  result  in  sa¬ 
tellites  in  the  photoelectron  spectra* ' *-6  and  in  the  Auger  and 
x-ray  spectra  from  transitions  through  which  the  photoexcited 
states  decay.2'4,8  The  study  of  these  multiple  excitation  proces¬ 
ses  is  important  because  they  epitomize  the  breakdown  of  the 
independent-particle  model  and  can  provide  important  clues  for 
the  understanding  of  electron  correlation  and  of  excitation  dyna¬ 
mics.  3 ' 5 ' 9 ' 10  The  energy  dependence  of  the  cross  sections  for 
double  excitation  is  particularly  informative  near  threshold;  the 
observation  of  Auger  satellites  makes  it  possible  to  measure  this 
dependence.  We  have  performed  an  investigation  in  which  highly 
monochromatized,  hard  synchrotron  radiation  was  tuned  through  the 
thresholds  for  various  multiple  excitation  processes  during  Is 
ionization  of  Ar,  and  the  probabilities  of  accompanying  3s  and  3p 
excitation  were  traced  by  measuring  the  intensities  of  pertinent 
Auger  satellites.  Results  were  compared  with  theory.3'*0 

In  the  experiment,  x  rays  from  an  8-pole  wiggler,  operating 
at  14  kG,  were  focused  onto  an  Ar  jet  by  a  Pt-coated  doubly 
curved  toroidal  mirror.  The  x-ray  bandwidth  from  a  Ge  (111) 
double-crystal  Bragg  monochromator  was  0.9  eV  at  hv=3200  eV;  the 
flux  on  target  was  -1012  photons/s  with  60  mA  of  3-GeV  electrons 
in  the  SPEAR  storage  ring.  Electron  spectra  were  measured  with  a 
computerized  double  cylindrical-mirror  analyzer;  with  a  pass 
energy  of  82.5  eV,  the  electron-spectrometer  resolution  was  1.6 
eV. 

We  take  the  2,660-eV  Ar  K-f,2L3  *D2  Auger-electron  lane  as 
reference.  The  K-LL  Auger  yield1-*  of  Ar  is  affected  only  minutely 
by  excitation  of  one  or  two  M-sheil  electrons.  The  intensity  of 
satellites  of  the  *D  Auger  line  relative  to  that  of  the  "diagram" 
line  is  therefore  a  measure  of  the  probaoility  of  the  multiple 
photoexcitation  processes  studied  in  this  work. 

To  interpret  the  *D  Auger  satellite  spectrum  it  was  necessa¬ 
ry  to  calculate  the  radiationless  transition  energies  and  rates 
in  the  presence  of  one  or  two  open  M  subshells.  The  initial 
states  can  be  limited  to  those  which  in  the  sudden  approximation 
are  expected  to  be  significantly  populated.1  These  are  the  [ls31] 
(1*0,1)  shakeoff  states  and  the  [ls31]nl  (n=4  and  5  for  1=1,  n=4 
for  1=0)  shakeup  states,  where  sauare  brackets  indicate  hole 
states.  In  the  limited  3s ( 2 S ) ns 1 ’ 3S ls2S  and  3p5 ( 2 P ) no1 ' 3S 1 s 2 S 
basis,  the  eigenstates  are  linear  superpositions  of  and  3S 
states.  The  initial  shakeup  states  can  be  identified  as  states 
with  dominant  *S  component  because  the  monopole  selection  rules 
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prevent  transitions  to  the  triplet  state.  According  to  our  Har- 
tree-Fock  (HF)  calculations,  the  initial  [ls3p]4p  state  has  al¬ 
most  pure  1S  character,  whereas  in  the  other  shakeup  cases  the 
mixture  is  more  uniform. 


FIG.  3.  Calculated  energies  of 
Auger  satellites  caused  by  3s-  and 
3p-electron  excitation  accompanying 
Is  ionization,  with  reference  to  an 
Ar  K-L2  3L2  3  Auger  spectrum  photoex- 
cited  2',  0  0d  eV  above  the  Is  ioniza¬ 
tion  threshold.  Estimates  of  relative 
satellite  intensities  within  each 
multiplet  are  indicated  by  the 
heights  of  the  bars. 


Calculated  Auger  satellite  energies  are  indicated  schemati¬ 
cally  in  Fig.  3.  The  satellites  arising  from  3s  and  3p  shakeoff 
accompanying  Is  ionization  are  seen  to  fall  into  the  peak  around 
-2,643  eV,  while  most  3s  and  3p  shakeup  processes  cause  Auger 
satellites  that  fall  within  the  2,650-eV  peak,  unresolved  from 
the  K-L 2I* 2  lso  diagram  line.  The  measured  intensity  of  the  *S 
line,  excited  below  the  threshold  for  any  [lsnl]  double  proces¬ 
ses,  is  11.0+0.6%  of  the  1D-line  intensity,  in  excellent  agree¬ 
ment  with  the  prediction  (11.12%)  from  a  relativistic  intermedia¬ 
te-coupling  calculation  that  includes  configuration  interac¬ 
tion.  The  predicted  positions  of  the  Auger  lines  are  only 
slightly  affected  by  configuration  mixing  in  the  initial  states 
and  by  relativity. 

In  Fig.  4(a),  the  relative  intensity  of  the  2,650-eV  shakeup 
satellite  peak  including  the  -^Sq  diacram  line  is  plotted.  The 
satellite  peak  that  arises  at  photon  ez ergy  E=3,225  eV  is  tenta¬ 
tively  ascribed  to  the  [ls3p]4p^  bound-bound  resonance,  in  accor¬ 
dance  with  the  interpretation  of  the  Ar  K  absorption  spectrum8 
which  shows  a  peak  at  3,224  eV.  In  accord  with  observations  of 
Kobrin  et  al.,6  the  [ls3p]4p  shakeup  satellite  appears  to  have 
approximately  half  its  asymptotic  intensity  at  5  eV  above  thres¬ 
hold.  Saturation  of  [ls3p]4p  plus  opening  of  the  [ls3s]4s  channel 
lead  to  a  small  gradual  increase  which  levels  off  to  a  constant 
shakeup  satellite  intensity  -60  eV  above  the  [ls3p]4p  threshold. 

If  we  assume  in  accordance  with  Fig.  4(a)  that  the  shakeup 
ratio  .is  practically  constant  as  a  function  of  E,  then  it  can  be 
concluded  that  the  experimental  shakeoff  curve  levels  off  at  high 
E  at  19+2%  (after  subtracting  the  threshold  value  of  5%).  The 
shape  of  the  curve  is  well  predicted  by  our  calculations  using  HF 
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wave  functions  for  both  the  [ls3p]  core  and  continuum  states 
[Fig.  4(b)].  In  the  calculation  of  the  continuum  wave  functions 
the  Lagrangian  multipliers  were  neglected,  but  a  Schmidt  orthogo- 
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FIG.  4.  (a)  Intensity  of 
the  2,650-eV  feature  in  the 
photoexcited  Ar  K-L2  3L2  3  Auger 
spectrum,  with  reference  to  the 
1D-line  intensity,  as  a  function 
of  x-ray  energy.  The  dashed  line 
at  11. IX  indicates  the  "S  diag¬ 
ram-line  contribution.  Energy 
thresholds  for  Is  ionization  and 
for  3  p  ->  4  d  and  3s*4s  shakeup 
accompanying  Is  ionization  are 
indicated  by  vertical  arrows. 
The  normalized  theoretical  pre¬ 
diction  for  the  near-threshold 
energy  dependence  of  these  rela¬ 
tive  shakeup  probabilities  is 
represented  by  the  solid  curve, 
(b)  Photoexcitation-energy  depe¬ 
ndence  of  the  2,643-eV  Auger 
satellite-group  intensity.  Thre¬ 
sholds  for  Is  ionization  alone 
and  accompanied  by  3p  and  3s 
ionization  are  indicated  by  ver¬ 
tical  arrows.  The  normalized 
theoretical  relative  shakeoff 
probability  is  indicated  by  the 
solid  curve.  Circles  and  trian¬ 
gles  pertain  to  data  from  sepa¬ 
rate  experiments. 


nalization  was  carried  out  afterwards.  As  seen  in  Fig.  4(b),  the 
measure}  cross-section  curve  is  only  slightly  affected  by  the 
opening  of  the  [ls3s]  shakeoff  channel.  Our  calculations  predict 
411  asymptotic  shakeoff  probability  of  25X  at  high  E. 

We  can  draw  the  following  conclusions.  (1)  The  difference  in 
the  photon-energy  dependence  of  shakeup  vs.  shakeoff  close  tc 
threshold  has  been  shown  experimentally,  for  the  first  time,  to 
be  as  predicted  by  theory.  (2)  The  measurements  indicate  more 
shakeoff  than  shakeup  at  high  photon  energy,  in  contrast  to  the 
predictions  of  Ref.  10  but  in  accord  with  Ref.  1.  (3)  The  mea¬ 
sured  shakeup  probabilities  agree  well  with  the  predictions  of 
Ref.  10,  but  the  shakeoff  probabilities  do  not.  (4)  The  measured 
total  shake  probabilities  are  bracketed  by  the  sudden-approxima¬ 
tion  values  calculated  by  the  HF  (DF)  method  for  a  double-hole 
[ls31]  and  a  single-hole  [Is]  field.  Within  the  restricted  HF 
(DF)  method,  both  procedures  are  somewhat  inconsistent  with  res¬ 
pect  to  fulfilling  the  closure  relation.  This  inconsistency  could 
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be  removed  if  many-e  1  ec tr on  wave  functions  were  used  to  obtain 
the 

<f(  [ls3  1]n(t)  D2SMffrozen([ls])2S> 

shakeup  and  shakeoff  amplitudes,  since  the  functions  are  eigen¬ 
functions  of  the  same  projected  (N-l) -electron  Hamiltonian  PH(N- 
1)P,  where  P*|ls><ls|. 

This  work  was  submitted  for  publication  during  the  period 
covered  by  the  present  Report  (cf.  Sec.  3,  Item  7)  and  has  since 
been  published  [Phvs.  Rev.  Lett.  54.  182  (1985)]. 
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2 . 3  Ionization  by  charged  particles 

Atomic  inner-shell  ionization  by  charged-particle  impact  has 
been  studied  intensively  during  recent  years.  The  theory  has 
progressed  beyond  the  plane-wave  Born  approximation,  through 
incorporation  of  binding  and  Coulomb-deflection  effects.1,2  For 
low-energy  collisions,  it  has  been  found  quite  important  to 
include  relativity  and  to  employ  more  realistic  wave  functions 
than  hydrogenic  ones.3-10 

Previously  tabulated  cross  sections  for  atomic  inner-shell 
ionization  by  charged-particle  impact  were  based  on  plane-wave 
Born-approximation  (PWBA)  calculations  with  hydrogenic  wave  fun¬ 
ctions.  1  12  No  systematic  calculations  of  inner-shell  ionization 
cross  sections  have  hitherto  been  published  that  take  into  ac- 
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count  the  effects  of  realistic  wave  functions  and  that  consisten¬ 
tly  incorporate  relativity. 

We  have  drawn  upon  our  previous  work  on  the  subject4,7'8  and 
invested  a  substantial  amount  of  computing  time  on  the  NASA  Ames 
Research  Center  CRAY  1  to  prepare  tables  of  relativistic  ioniza¬ 
tion  cross  sections  for  iC  and  L  shel  Is  of  27  elements  with  atomic 
numbers  22<Z<92,  for  proton  impact  with  incident  energies  from 
0.15  to  3  MeV.  These  cross  sections  were  calculated  with  Dirac- 
Hartree-Slater  wave  functions  and  include  binding,  polarization, 
and  Coulomb-deflection  corrections. 

The  work  has  been  accepted  for  publication  in  Atomic  Data 
and  Nuclear  Data  Tables  (cf.  Sec.  3,  Item  10). 
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2.  W.  Brandt  and  G.  Lapicki ,  Phys .  Rev.  A  23,  1717  (1981). 

3.  G.  Lapicki  and  A.  R.  Zander,  Phys.  Rev.  A  23,  2072  (1981). 

4.  M.  H.  Chen.  B.  Crasemann,  and  H.  Mark,  Phys.  Rev.  A  26,  1  243 
(1982). 

5.  T.  Mukoyama  and  L.  Sarkadi ,  Phys.  Rev.  A  23,  375  (1981). 

6.  T.  Mukoyama  and  L. Sarkadi,  Phys.  Rev.  A  25.  1411  (1982). 

7.  M.  H.  Chen,  B.  Crasemann,  and  H.  Mark,  Phys.  Rev.  A  27,  2358 
(1983). 

8.  M.  H.  Chen,  Phys.  Rev.  A  30,  2082  (1984). 

9.  0.  Benka  and  A.  Kropf,  At.  Data  Nucl.  Data  Tables  22,  219 

( 1978) . 
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2.4  Lifetimes  of  deep  hole  states 

Atomic  inner-shell  vacancy  states  decay  predominantly  by 
radiationless  processes  and,  to  a  minor  extent,  by  radiative 
transitions.  The  coupling  between  radiative  and  radiationless 
channels  has  usually  been  neglected.1  The  electrostatic  coupling 
between  the  final  K-LL  Auger  channels  in  Ne  was  considered  by 
Howat  et  al.,2'3  and  for  K-LL,  K-LM,  and  K-MM  transitions  in  Mg 
by  Howat.4  These  authors  found  that  interchannel  coupling  substa¬ 
ntially  alters  the  spectral  distribution  of  the  K  Auger  electrons 
which  are  emitted  with  kinetic  energies  of  several  hundred  elec¬ 
tron  volts.  It  is  tempting,  therefore,  to  investigate  similar 
effects  in  the  case  of  Coster-Kronig  processes  which  are  extreme¬ 
ly  intense,  compared  with  Auger  processes,  and  in  which  the 
outgoing  electrons'  kinetic  energy  is  at  least  an  order  of  magni¬ 
tude  less  than  in  Auger  transitions. 
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One  inviting  case  for  theoretical  investigation  is  that  of 
L-shell  Coster-Kronig  transitions  in  Ar,  for  which  perplexing 
discrepancies  have  existed  between  experimental  results5,6  and 
theoretical  Dredictions.7*11  Hartree.'  Hartree-Fock-S later  ,8  and 
Hartree-Fock9- 10  single-configuration  calculations  have  been 
found  to  overestimate  experimental  A.-  L.-L23M.  Coster-Kronig 
rates  by  factors  of  -4,  the  L1“L23M1  to  L1-L,^M,3  intensity  ratio 
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Possible  reasons  for  these  striking  discrepancies  in  theore¬ 
tical  predictions  include  (i)  many-body  interactions  in  the  ini¬ 
tial  and  final  atomic  systems,  (ii)  effects  of  relaxation  in  the 
final  ionic  state,  and  (iii)  effects  of  the  exchange  interaction 
between  the  continuum  electron  and  the  final  bound-state  elec¬ 
trons.  These  possibilities  have  been  investigated  in  recent  cal¬ 
culations  . 


Mehlhorn1  has  recently  reevaluated  his  experimental  data, 
leading  to  a  13*  decrease  in  the  L x -L2 3 M 1 / M 2 3  intensity  ratio. 
This  revised  experimental  ratio  agrees  better  with  theoretical 
predictions,  but  the  discrepancies  remain  large  in  comparison 
with  other  atomic  transition-rate  calculations  in  which  theory 
tends  to  come  within  10-20*  of  measured  rates. 

In  an  effort  to  resolve  this  perplexing  situation,  we  have 
investigated  the  effects  of  continuum  interaction  in  the  outgoing 
channels  on  the  total  and  relative  term  intensities  of  Ar  L-shell 
Coster-Kronig  transitions.  The  rates  were  calculated  through  a 
formulation  which  includes  the  effects  of  final-state  channel 
mixing.  The  interchannel  interaction  was  found  to  change  indivi¬ 
dual  transition  rates  by  as  much  as  84*.  The  relative  term  inten¬ 
sities  from  this  calculation  agree  better  with  experiment  than 
computations  which  do  not  take  account  of  channel  mixing.  The 
discrepancies  previously  found  between  calculated  and  measured 
total  Coster-Kronig  rates  and  between  1P/3?  ratios  are  largely 
removed  when  channel  coupling  is  taken  into  consideration. 

A  description  of  this  research  has  been  accepted  for 
publication  in  Physical  Review  A  (cf.  Sec.  3,  Item  11);  earlier 
work  on  this  subject  has  been  published  in  the  course  of  the 
report  year  (Sec.  3,  Items  2,  4,  6). 
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2 . 5  X-ray  emission  from  the  M  shell 

Most  of  the  existing  calculations  of  x-ray  emission  rates 
deal  with  transitions  to  K-  and  L-shell  vacancies.  Only  one 
calculation  has  been  performed  for  M  x-ray  emission  rates,  based 
on  relativistic  Har tree-S 1 ater  wave  functions  and  covering  6 
elements  with  atomic  numbers  between  48  and  93.  We  have  now 
performed  Dirac-Fock  calculations  of  M  x-ray  emission  rates,  both 
in  the  Coulomb  and  length  gauges,  for  10  elements  with  atomic 
numbers  between  48  and  92.  We  have  studied  the  gauge  dependence 
of  the  M  x-ray  rates,  which  is  due  to  the  effect  of  the  nonlocal 
potential.  We  have  explored  the  effects  of  relativity  on  these 
rates,  which  are  caused  by  changes  in  transition  energies,  shifts 
in  the  wave  functions,  and  the  role  played  by  higher  multipoles. 
Typical  results  are  illustrated  in  Fig.  5.  This  research  has  been 
published  (Section  3,  Item  3.) 
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FIG.  5.  Total  probabilities  (in 
atomic  units)  of  radiative  transi¬ 
tions  to  an  M ^  vacancy,  as  functions 
of  atomic  number  Z.  Results  from 
relativistic  Dirac-Fock  calculations 
in  the  length  gauge  (DF-L)  and  in  the 
Coulomb  gauge  (DF-C)  are  compared 
with  rates  from  relativistic  local- 
potential  Dirac-Hartree-Slater  (DHS) 
and  nonrelativistic  Har t ree-S 1  at e r 
(HS)  calculations . 
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